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This paper lays out a comprehensive methodology for computing a low-speed, high-lift 


available at the conceptual design stage. Introducing low-order, physics-based aerodynamic 
analyses allows the methodology to be more applicable to unconventional aircraft concepts 
than traditional, fully-empirical methods. The methodology uses empirical relationships 
for flap lift effectiveness, chord extension, drag-coefficient increment and maximum lift 
coefficient of various types of flap systems as a function of flap deflection, and combines 
these increments with the characteristics of the unflapped airfoils. Once the aerodynamic 
characteristics of the flapped sections are known, a vortex-lattice analysis calculates the 
three-dimensional lift, drag and moment coefficients of the whole aircraft configuration. 
This paper details the results of two validation cases: a supercritical airfoil model with 
several types of flaps; and a 12-foot, full-span aircraft model with slats and double-slotted 
flaps. 


T akeoff and landing performance requirements often place strong constraints on sizing of the wing and 
engines during conceptual aircraft design, yet the low-speed aerodynamic characteristics of an aircraft 
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polar, without requiring additional details about the aircraft design beyond what is typically 
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I. Introduction 
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with its flaps deflected and landing gear extended can be difficult to estimate accurately. In particular, 
the maximum lift coefficient is a difficult value to obtain accurately, even when comprehensive information 
about the flap and slat geometry is available and the configuration is analyzed using higher-order analysis 
methods.^’ ^ Traditionally, designers have tended to rely on scaling of available historical data and test 
results,^ but these are of limited use with newer unconventional configurations for which no full-scale data 
are available. The research described herein seeks to fill this capability gap by providing a comprehensive 
method for computing a low-speed, high-lift polar without greatly expanding the amount of detail about 
the aircraft design required to perform the analysis. Introducing low-order, physics-based aerodynamic 
analyses allows the methodology to be more applicable to unconventional aircraft concepts than traditional, 
fully-empirical methods. 


II. Methodology 

The methodology for the aerodynamics of an aircraft in the high-lift configuration involves building up the 
total lift and drag curves in a quasi-three-dimensional manner. First, the shape of the two-dimensional section 
lift-curve slope is determined from base (unflapped) wing sectional characteristics; second, two-dimensional 
increments to the section lift-curve slope due to flap deflection are estimated based on information about 
the flap type and basic geometry; and third, the three-dimensional aerodynamic coefficients of the complete 
aircraft are determined based on the base and incremental aerodynamic characteristics of the wing sections. 
The following sections describe these stages in detail. 

A. Base Section Aerodynamics 

Figure 1 shows an idealized (linear) section lift curve; in reality the slope of the curve would decrease at 
higher angles of attack approaching the maximum lift coefficient. However, even though takeoff and landing 
performance is strongly constrained by the maximum lift coefficient of the aircraft, the actual low-speed 
operation of the aircraft takes place at more benign lift coefficients where the lift curve is nearly linear. 
Using this assumption, one can represent the lift curve of the section by three parameters: the lift-curve 
slope, the angle of attack at zero lift, <a/=o (or alternatively, the zero-angle lift coefficient, Qq), and 
the maximum lift coefficient, Qj^iax* DATCOM^ provides separate semi-empirical methods for the 

calculation of these parameters. 



Figure 1. Idealized section lift curve 

Alternatively, the section aerodynamics may be determined directly using computational methods. In 
this study, two different computational methods were evaluated for this purpose: XFoil,^ a two-dimensional 
panel method for single-element airfoils coupled with an integral boundary-layer (IBL) analysis, and MSES,^ 
a two-dimensional Euler analysis for single- and multiple-element airfoils also coupled with an IBL analysis. 
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Section III provides a comparison of DATCOM versus computational methods in the analysis of sectional 
characteristics. 

B. Increments Due to Flap Deflections 

When a trailing-edge flap is extended, the following changes occur to the section lift-curve slope, as shown 
in Fig. 2: the lift at zero angle of attack is increased by an increment AjQq, the lift-curve slope changes 
to a modified value of and the maximum sectional lift coefficient increases by an increment 
Methods for estimation of these three effects will be investigated in the following sections. Except where 
otherwise noted, the methodology in this section follows the approach of Ref. 7. 



Figure 2. Two-dimensional trailing-edge flap effects 


1. Lift Increment at Zero Angle of Attack 

Flap lift effectiveness is characterized as the change in lift coefficient at zero angle of attack due to flap 
deflection. Glauert’s linearized theory for thin airfoils with flaps^ gives the following result for the theoretical 
increment due to a given trailing-edge flap deflection, Sf: 


(AjQo)theory 

where as is the rate of change of zero-lift angle of attack with flap deflection, defined as follows: 



( 1 ) 

( 2 ) 


where ci^ = ^ is the rate of change of lift with flap deflection at constant angle of attack and ^ is 

the rate of change of lift with angle of attack at constant flap deflection. For a thin airfoil, as is related to 
the flap geometry as follows: 


Of — sin Of 
as = 1- ^ 

7T 


(3) 


where 

= cos-1 


(4) 
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Here c/ is the flap chord. 

In practice, the theoretical lift increment in Eq. (1) cannot be realized because of several non-linear 
effects. First, the small-angle assumption of linear theory breaks down at large flap deflections. Second, the 
de-cambering effect of the boundary layer near the trailing edge of the flap results in a smaller lift coefficient 
than predicted.^ Finally, at large flap angles the flow can become separated. These effects are captured 
using an empirical flap effectiveness factor, rjs, which is applied to the theoretical equation: 

A/Cio = ijgascijf (5) 

Figure 3 shows empirical curves of flap effectiveness factor for various types of flap systems as a function 
of deflection angle. These curves allow the conceptual designer to calculate the lift characteristics of repre- 
sentative flap systems under real-world conditions, but they do not require detailed information about the 
geometry of the flaps. 

In addition to changing the camber and incidence of the base airfoil, slotted and Fowler flaps also have 
the benefit of increasing the effective chord of the section through rearward motion of the flap segments. 
Since the circulation has a larger chord length to operate on, the flap extension results in an additional 
increase in lift coefficient proportional to the ratio of the new and old chords. If one assumes that the base 
section ci^ is not altered by the chord extension, the lift coefficient increment with flaps extended can be 
determined for the extended chord and then referenced against the original chord as follows: 

^/Cio = ^fc'lo 7 + (7 “ 

where c' = c -h Ajc is the extended chord and is the lift increment referenced against the extended 

chord: 

^f^io = '^sa'sCi^Sf (7) 

The value of a'^ is found from Eqs. (3) and (4) by replacing the original flap chord with the extended flap 
chord, Cj = cj -h Ajc. If known, the extended chord can be used directly in Eq. (6), or may be estimated 
from the empirical curves shown in Fig. 4. 

2. Lift- Curve Slope 

Deflection of trailing-edge flaps has a number of different effects on the lift-curve slope of a section. The 
potential-flow effect of flap deflection decreases with increasing angle of attack, while viscous effects cause a 
decrease in flap effectiveness at larger angles of attack. Finally, the chord extension of slotted flaps causes 
a proportional increase in q^. These effects are complex and difficult to determine accurately; however, a 
rough estimate of the effects is provided in Ref. 7: 



This equation is appropriate for smaller angles of attack (0 < a < 5°), provided that the slot shape is well 
optimized and stall has not occurred. 

3. Maximum Lift Coeffieient 

The effect of trailing-edge flap deflection on wing maximum sectional lift coefficient is strongly influenced by 
the stalling behavior of the base airfoil section. This stalling behavior is closely correlated with the leading- 
edge sharpness parameter, A^, which is the difference between the upper-surface ordinates at 6% chord and 
0.15% chord, respectively, normalized by chord. Airfoils with sharp noses (A^ up to 1.5% chord) exhibit 
leading-edge, or long-hubble stall behavior, in which a region of separated flow at the leading edge progresses 
rearward as the angle of attack is increased, leading to stalling of the entire section. For these sections, 
thin-airfoil theory finds that the increment in maximum lift coefficient due to flap deflection is related to the 
theoretical lift increment: 

/^^max = ( /^^oJtheory _gin(9/) ^ ^ 
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(a) Plain flaps with sealed gap 




Figure 3. Empirical lift effectiveness, rjs for various types of trailing-edge flaps as a function of deflection angle. 
Source: data from Ref. 7. 
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Figure 4. 


Flap chord extension as a function of deflection for different flap types. Source: data from Ref. 7. 


For practical ranges of flap chord ratio, Eq. (9) results in a which is about half of the theoretical lift 

increment. Combining Eq. (9) with Eqs. (4) and (7) yields the following simple expression for the maximum 
lift increment for airfoils with sharp noses {Ay < 1.5): 


^•f^^max 


C-lJf 


sin Of 

7T 


( 10 ) 


As Ay increases above 1.5, the section displays a short-bubble stalling behavior and the maximum lift 
increment increases approximately proportionally with Ay to the point where the stall is dominated by 
trailing-edge stall. Beyond this value of A^, the base section displays a separation behavior that is associated 
with the pressure gradient near the trailing edge; therefore, the stall will be delayed by the local suction 
produced by the initial deflection of a trailing-edge flap.^^ Beyond a certain flap angle, however, the load 
induced at the airfoil nose will be high enough that stall will occur first at the leading edge and the maximum 
lift increment will behave in a manner similar to the sharp- nose airfoil relationship in Eq. (10). 

The maximum lift coefficient for airfoils with short-bubble or trailing-edge stalling behavior {Ay > 1.5) 
may thus be determined by finding which of the two stalling mechanisms is dominant: 


c 


toax 



6 f=o ^f^io 

/ X 0.08 

(sxioe) 2 (^^0 + ^f^lo) 


( 11 ) 


The top term in Eq. (11) reflects the one-to-one relationship between the lift increment and the increase in 
maximum lift for airfoils experiencing trailing-edge stall at smaller flap deflections, whereas the lower term 
reflects the increase in maximum lift for airfoils experiencing leading-edge stall at larger flap deflections. 
Ref. 7 suggests that this equation can be expected to give acceptable results for both plain and slotted flaps, 
including multi-element flaps with or without Eowler motion. 


4 . Profile- Drag Inerement 


Basing AfCi on Glauert’s linear airfoil theory and accounting for chord extension, the minimum-profile-drag 
increment can be expressed as 


^f^dpo 


C f 

kdCi„a's—5f sin (5/ 


■ Qp 



( 12 ) 


where a'^ is defined as before. The profile-drag factor kd is an empirical function of the type of flap and the 
deflection angle, as shown in Eig. 5. This figure shows that the actual drag increment is on the order of 20 
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Figure 5. 


Effect of trailing-edge flap deflection on profile drag factor, kd- Source: data from Ref. 7. 


to 25% of the theoretical drag increment for plain flaps, and 8 to 15% for slotted flaps. The high value of kd 
for slotted flaps at small deflection angles is due to unfavorable slot flow and depends to a large extent on 
the details of the flap system design. 


5. Effects of Leading-Edge Elaps 


Even airfoils with rounded noses will exhibit leading-edge stall at larger trailing-edge flap deflection angles. 
When leading-edge devices (either flaps or slats) are deployed, the high suction peak near the leading edge 
of the airfoil is alleviated, delaying the leading-edge separation to a higher angle of attack and increasing 
the maximum lift coefficient.^^ However, a slat wake may interfere with the flow around trailing-edge flaps 
and the increment in Qj^iax than for the undeflected case. 

Leading-edge droop causes a loss in lift at zero angle of attack which can be derived from Glauert’s thin 
airfoil theory: 


^sClo 


0a - sin 0a s- 

^S^loc 

7T 


(13) 


where 

Os = cos“^ — 2—^ (14) 

If one assumes that the loss of lift is entirely due to a change in the ideal angle of attack of the section — that 
is, the angle at which the flow approaches the leading edge of the airfoil smoothly and causes no extreme 
suction peak at the nose — then the change in maximum lift coefficient can be derived from DATCOM data 
to give 




max 



(15) 


This relationship is valid for deflection angles up to approximately 25 degrees. Leading-edge flap deflection 
has no appreciable effect on the lift-curve slope provided that the flap is not deflected much more than 20 
degrees. 

As in the case of plain leading-edge flaps, extension of a slat results in a loss of lift at zero angle of attack. 
However, a slat also increases the effective chord length; both effects are of the same order of magnitude and 
one may assume that they cancel each other out. Assuming that the trailing-edge flap and the slat have been 
designed together so that all parts of the system experience separation at nearly the same angle of attack, 
the following expression may be used for the airfoil with slats deployed: 


^^max 


(1 - ka) 


Qo + ^/Qo + 0-47c;^ 
1 + 0.035c;^ 


( 16 ) 
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where kg is a factor taking into account the non-linearity of the lift curve at high lift and generally varies 
between 0.03 and 0.15, with 0.07 as a reasonable average. The maximum lift of Krueger flaps is very similar 
to that of slats and Eq. (16) may be used in this case as well. 

C. Three-Dimensional Aerodynamics 

The final step in the low-speed analysis is to calculate the three-dimensional aerodynamic coefficients of the 
complete aircraft based on the base and incremental aerodynamic characteristics of the wing sections. In 
this methodology the Athena Vortex-Lattice (AVL) computer program^ ^ is used. Given the two-dimensional 
characteristics of the wing and tail sections, AVL can be run with either an ideal lift curve slope = 27ra) 
or with a lift-curve slope determined from DATCOM or from two-dimensional computational analysis; the 
effect of this choice is examined in Section III. 

AVL models all lifting surfaces as “flat”; all sections of a surface are parallel to the x-axis, and the 
incidence and camber of lifting surfaces are modeled as a tilting of the control point normal vectors rather 
than being reflected in the physical shape of the surface. Along these lines, the deflection of control surfaces 
(such as flaps, elevators, and ailerons) is also modeled as a tilting of control point normal vectors rather than 
a physical deflection of the surface. AVL is a linear code, so flap deflections produce lift increments in a 
manner consistent with Eq. (1). To account for the non-linear effects discussed previously, the empirical lift 
effectiveness factor, r]s from Eq. (5), is determined for each flap component and input to AVL as the ratio 
between the effective and actual flap deflections, referred to in the program as flap gain. The increment in 
zero-lift drag is also added to the section information, while the increment in pitching-moment is calculated 
by AVL. 

The chordwise extension of slotted flaps is modeled as a physical chord increase in the AVL geometry, 
but for calculation of aerodynamic coefficients the reference area remains unaffected. This chord extension 
can come from either historical data for different flaps, as seen in Eig. 4, or from the actual chord extension 
of the flap, if known. Since in this methodology slats have no influence on lift, they are not modeled in 
the AVL geometry, but are factored into the zero-lift drag of the section and are used in the calculation of 
maximum sectional lift coefficient. 

In addition to the lift-curve slope, the zero-lift drag coefficient from DATCOM and the flap drag increment 
from Eq. (12) are input to AVL for each section of the wing and tails, with the program automatically 
integrating the values to determine the total profile drag of each lifting surface. The methodology calculates 
the drag of the fuselage using the ERICTION computer program^ ^ — which determines the profile drag of 
aircraft components using empirical skin-friction drag of a fiat plate combined with empirical form factors — 
and adds it to the drag of the lifting surfaces to arrive at the drag of the complete aircraft. 

The maximum lift coefficient of the aircraft, calculated using the critical section methodd^ 

for increasing values of angle of attack, the lift coefficient at each location along the wing is compared to 
the maximum lift coefficient of that section and a stall is declared when the maximum lift coefficient is first 
surpassed at any section. The maximum lift coefficient for the aircraft is defined as the value achieved when 
the stalling condition is reached. The critical section method is simple and computationally inexpensive; for 
a highly-swept wing, however, the actual ^^^ds to be lower due to the effects of spanwise flow at 

higher angles of attack, and studies have found that the critical section method can be optimistic in these 
cases. 


III. Validation 

Validation of the methodology was performed for two different cases. Eirst, a supercritical airfoil model 
was used to test the validity of the two-dimensional aerodynamic analysis of the base section and associated 
increments due to flap and slat deflection. Second, a full-span, high-aspect-ratio transport model was used 
to test the validity of the three-dimensional aerodynamic analysis. 

A. Supercritical Airfoil with Flaps 

The two-dimensional analysis methodology was validated using experimental results for a 9.3%-thick su- 
percritical airfoil. Test results are available for this airfoil for four different geometric configurations: 
the base airfoil without flaps or slats installed, and the same airfoil configured with a single-slotted flap, a 
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double- slotted flap and slat combination, and a triple-slotted flap and slat combination. In this study the 
first three configurations were analyzed. 

Figure 6(a) shows analysis results for the base airfoil using three different methods: DATCOM empirical 
predictions, and MSES and XFoil analyses. The DATCOM predictions show good agreement for the lift- 
curve slope of the airfoil but under-predict the zero- angle lift coefficient and the maximum lift coefficient. 
The DATCOM database is based on earlier airfoil designs, and newer supercritical airfoils tend to fall outside 
the range of some of the correlating parameters, due in particular to the large amount of aft camber. The 
estimate of zero-lift drag compares favorably with the experimental data, but DATCOM does not provide 
methods for profile drag-due-to-lift or moment coefficient so those comparisons are not included. 

Both MSES and XEoil give accurate estimates of the lift curve, drag polar and moment-coefficient curves 
for the base airfoil. XEoil, being a panel code, is robust at high angles of attack and gives results up to and 
past stall, whereas MSES tends to only return results up to the vicinity of the stall. Note that the data 
points in Eig. 6(a) do not show a clear stalling point, but the final data point exhibits a rapid drag increase 
and appears to be close to the stall. If this is true, XEoil estimates maximum lift coefficient accurately for 
this case, whereas MSES slightly under-predict s it. 

Analysis of the airfoil with a single-slotted flap was performed using XEoil results for the base airfoil 
combined with the empirical flap increments described previously, as well as using direct analysis of the 
airfoil and flap geometry in MSES (Eig. 6(b)). Both the empirical curves and MSES give good results for the 
lift increment due to flap deflection. Both methods predict a maximum lift coefficient just over 2.5 but again 
the data do not indicate a clear stalling point so the experimental maximum lift coefficient is unknown. The 
drag coefficient increment added to the XFoil analysis gives a reasonable estimate of the total drag with the 
flap deflected, but since it is not a function of angle of attack the shape of the drag polar differs from the 
experimental results. Both analysis methods over-estimate the nose-down pitching- moment increment of the 
trailing-edge flap deflection. The ability of the two-dimensional methods to calculate the moment-coefficient 
increment, however, is not relevant here because this increment is calculated by AVL in the full methodology. 

Finally, the airfoil with double-slotted flap and slat was analyzed using XFoil only, combined with the 
empirical flap increments, since MSES would not converge for this case (Fig. 6(c)). The empirical curves give 
an accurate estimate of the lift increment, but somewhat overestimate the drag increment due to flap deflec- 
tion. Contrary to the single-slotted flap case, the moment-coefficient increment is greatly underestimated 
for the double-slotted case. 

Table 1 gives a summary of the maximum lift coefficient calculated by the different methods compared 
to the experimental results for each case. The “> ” symbol indicates values for which the actual stall point 
is not clear from the data. The combination of XFoil with empirical flap increments provides reasonably 
accurate results for cases where the maximum lift coefficient can be determined from the data. 

Table 1. Maximum lift coefficient results — 9.3% supercritical airfoil, Mach 0.201, Rcc = 2.83 x 10® 


Analysis method 

Base airfoil 

Single-slotted 

Double-slotted w/ slat 

DATCOM 

1.53 

— 

— 

MSES 

> 1.61 

2.57 

— 

XFoil 

1.74 

2.51 

4.49 

Data 

1.72 

> 2.42 

> 4.57 


B. EET AR12 Model 

As a validation case for the three-dimensional analysis of an aircraft with flaps, the EET AR12^^’^^ model 
was used. This model is a 12-foot full-span aircraft configuration with a supercritical wing, full-span slats, 
and part-span double-slotted flaps with a cutout for the engine (Fig. 7). The model also can be equipped 
with a movable horizontal tail, flow-through nacelles, and landing gear, but these were not installed for the 
test cases used in this study. 

An OpenVSP^^ model of the EET AR12 was constructed using the published planform shape, airfoil 
coordinates, twist distribution, and fuselage cross-section shapes (Fig. 8). The OpenVSP model was then 
converted to an AVL model using an in-house automated conversion routine; flap and slat definitions were 
added to the AVL model during this conversion process. Since both the chord extension and flap gain are 
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Figure 6. Analysis results — 9.3% supercritical airfoil, Mach 0.201, Rcc = 2.83 x 10® 
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Figure 7. EET AR12 model planform detail. Source: Ref. 18 
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o.nflri 



Figure 8. EET AR12 OpenVSP model 


functions of the trailing-edge deflection, separate AVL models were required for the cruise wing (flaps stowed) 
and the takeoff and landing configurations, as illustrated in Fig. 9. The AVL model uses a cruciform fuselage 
model with no wake shedding that produces a moment in both the longitudinal and directional axes but 
near-zero lift and side forces for the fuselage. The wing surface includes a connecting section that carries 
the lift from the side of body across the fuselage centerline but is uncambered to correct for the poor lifting 
efficiency of the fuselage. The combination of a non-lifting fuselage with a wing carry-through was found 
to produce the best match between the predicted and experimental lift curves compared to other modeling 
strategies. 

Validation results are shown in Fig. 10 for the three cases: cruise, takeoff and landing. In each case the 
incidence angle of the horizontal tail was kept fixed at zero degrees. Since the vortex- lattice analysis is linear, 
the lift curves do not capture the effects of stall. However, the maximum lift coefficient of the configuration 
was calculated using the critical section method, with the overall lift curve calculated by the combination of 
DATCOM plus AVL, and with the sectional maximum lift coefficients calculated by XFoil. Table 2 compares 
the results of this analysis with experimental measurements. 


Table 2. Maximum lift coefficient results — EET AR12 model, Mach 0.168, Rcc = 137 x 10® 


Analysis method 

Cruise wing 

Takeoff flaps 

Landing flaps 

AVL + DATCOM + XFoil 

1.79 

2.41 

2.77 

Data 

1.34 

2.51 

2.82 


For the cruise wing (Fig. 10(a)), the lift curve and drag polar match well in the pre-stall region. The 
pitching- moment prediction is not as good and there is significant difference in the slope of the curve compared 
to experimental results. The maximum lift coefficient is over-predicted by 33%, but as noted previously the 
critical section method can be expected to be optimistic for a swept wing. 

For the takeoff configuration (Fig. 10(b)), the predicted lift curve matches well at moderate angles of 
attack, but there is a much larger difference at lower angles of attack. In the experimental data there is an 
unexpectedly low lift increment and non-linearity in the lift curve at near-zero angles of attack, suggesting 
that the combination of takeoff flap geometry and deflection was not particularly ideal at those conditions. 
The differences between the predicted and experimental lift curves in this region can therefore be expected, 
because the empirical curves for flap effectiveness assume a well-optimized flap design. The drag polar un- 
derestimates the induced-drag effects of the flap deflection, but the method gives a reasonable approximation 
of the drag coefficient at moderate angles of attack. The prediction of maximum lift coefficient compares 
favorably with the data. At moderate angles of attack, the slope of the pitching-moment curve is inaccurate 
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(a) Cruise wing 



(b) Takeoff wing 



(c) Landing wing 


Figure 9. EET AR12 AVL models for different flap settings 
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Figure 10. Analysis results — EET AR12 model, Mach 0.168, Rcc = 137 x 10® 
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despite the lift and drag being well-predicted. 

The landing configuration comparisons are shown in Fig. 10(c). Again, the lift curve and drag polar 
show good comparison with experimental results, and the prediction of maximum lift coefficient compares 
favorably with the experimental value. As before, the magnitude of the slope of the pitching-moment curve 
is significantly over-predicted relative to the data. 

IV. Conclusion 

The classical DATCOM predictions proved to be somewhat inadequate for the prediction of the lift curve 
and maximum lift coefficient of the airfoils used in this study. The predictions are based on empirical curve 
fits to airfoil existing at the time, and newer supercritical airfoil designs are not part of the database. Shape 
parameters for supercritical airfoils, particularly the location of the maximum camber, are not properly 
captured. Analysis of sectional aerodynamics using XFoil produced good agreement with experimental data; 
being a panel method, it is very robust at high angles of attack and predicts fairly well with the 

caveat that the onset of stall tends to occur more gradually than in the experimental data. Analysis with 
MSES also produced good agreement with experimental data but it was more difficult to produce converged 
solutions at high angle of attack in order to capture the stall point. 

For prediction of the aerodynamics of the full aircraft in the cruise configuration, vortex- lattice analysis 
compared favorably with experimental data. The lift curve was accurately modeled although the predicted 
C^max optimistic; this was as expected due to the importance of spanwise flow near stall for a swept 
wing. Both the zero-lift drag coefficient and the induced drag were matched well. The moment coefficient 
of the full configuration at zero lift was fairly accurate but the slope of the moment coefficient curve with 
angle of attack was not well predicted. 

For the takeoff flap setting, the increment due to flap deflection was predicted well except at low angles 
of attack where the experimental results seem to indicate a non-linearity that is not captured in the cur- 
rent methodology. The induced-drag effects of flap deflection were somewhat under-predicted; in contrast, 
moment-coefficient effects were notably over-predicted. For the landing flap setting, the predictions of the 
lift curve and drag polar compared well to the data, although the moment-coefficient increments were again 
over-predicted. For both the takeoff and landing flap settings, the predicted C^Lj^ax compared well with the 
experimental data. 

Plans exist to extend this methodology using higher-order aerodynamic methods such as a panel method 
and Euler analysis. Similar to the current methodology, empirical flap effectiveness factors can be used to 
define an effective flap deflection that takes into account the viscous effects without the need for a fully- 
viscous solution. The advantage of a full three-dimensional analysis is that it would enable the use of the 
pressure-difference rule for the prediction of C^Lj^iax’ which has been shown to produce robust and accurate 
predictions of wing stall. The primary disadvantage of a three-dimensional analysis would be that the flap 
and slat surface shapes would need to be defined in detail, rather than the more basic information, such as 
flap chord fraction, that is needed to perform a vortex-lattice analysis. 
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